clc,clear
close all
EMao = 30e-3;
F = zeros(11,3);
IndexF = 1;
StepDeltaX = EMao/10;
for DeltaX = 0:EMao/10:EMao
    F(IndexF,:) = FCase12(DeltaX);
    IndexF = IndexF+1;
end
subplot(1,3,1);
plot([0:EMao/10:EMao],F(:,1));
xlabel("Case1DeltaX");
ylabel("Case1Fx");
title("Fx-x");
subplot(1,3,2);
plot([0:EMao/10:EMao],F(:,2));
xlabel("Case1DeltaX");
ylabel("Case1Fy");
title("Fy-x");
subplot(1,3,3);
plot([0:EMao/10:EMao],F(:,3));
ylabel("Case1Fz");
xlabel("Case1DeltaX");
title("Fz-x");
